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UNCLASSIEIED 


Coulomb  Crystals  in  a 
Pulse-Excited  Linear  Paul  Trap 

N.  Kjajrgaard,  K.  M0lhave,  and  M.  Drewsen 


Institute  of  Physics  and  Astronomy,  Aarhus  University,  DK-8000  Aarhus  C,  Denmark 

Abstract.  We  present  results  from  recent  experiments  with  laser  cooled  ions  captured  in 

a pulse-excited  linear  Paul  trap,  As  discussed  in  [l-.'l],  the  radial/transverse  dynamics  of  charged 
particles  in  an  electric  quadrupole  trap  with  a pulsed  voltage  excitation  is  analogous  to  that  of 
particles  travelling  through  the  alternating  gradient  (AG)  magnetic  quadrupole  lattice  of  a storage 
ring  or  transport  channel.  Hence,  the  pulse-excited  Paul  quadrupole  trap  may  be  used  to  emulate 
non-neutral  beam  propagation.  In  a series  of  experiments,  we  have  investigated  the  stability  of 
stationary  Coulomb  crystals.  The  results  comply  with  current  theory  of  crystalline  plasmas  and 
beams.  We  observe,  however,  that  the  plasma  may  stay  crystalline  beyond  the  expected  stability 
region  in  a hitherto  unseen  configuration  consistent  with  a string  of  rotating  discs. 


INTRODUCTION 

When  a charged  particle  travels  through  a magnetic  quadrupole  it  experiences  a focusing 
force  in  one  dimension  and  a defocusing  force  in  the  perpendicular  dimension.  However, 
by  combining  two  quadrupoles  with  opposite  polarities  into  a doublet  a net  focusing 
effect  can  be  obtained  in  both  dimensions  as  illustrated  in  Fig.  1 . Alternating  gradient 
(AG)  magnetic  quadrupole  focusing  was  proposed  by  Courant,  Livingston,  and  Snyder 
[4]  in  the  early  1950’s‘  and  despite  its  simplicity,  this  concept  revolutionized  accelerator 
design  [6]  and  is  today  widely  used  in  storage  rings  to  achieve  transverse  confinement 
of  particles  along  the  design  orbit.  The  resulting  transverse  dynamics  can  be  emulated  in 
the  stationary  frame  of  en  electric  quadrupole  trap  [1-3]  and  in  recent  experiments  we 
have  exploited  this  analogy  to  investigate  the  stability  of  laser  cooled  ^'’Mg^  Coulomb 
crystals  produced  and  imaged  as  described  in  [7]. 


THE  PULSE-EXCITED  LINEAR  PAUL  TRAP 

The  linear  Paul  trap  employed  in  the  experiments  consists  of  four  cylindrical  electrodes, 
as  shown  in  Fig.  2(a),  with  a radius  of  /?  = 4.0  mm;  the  radius  of  the  inscribed  inter- 
electrode space  is  rq  = 3.5  mm.  To  operate  the  trap  in  a pulsed  mode,  the  periodic 
rf  voltage  (|).t(0  shown  in  Fig.  2(b)  is  applied  to  the  rods.  The  pulsed  waveform  is 
characterized  by  an  amplitude  Vrf,  a period  T,  and  a pulse  duration  xT.  As  indicated 


' The  idea  had  been  suggested  earlier  by  Christofilos  [5]  but  was  unpublished. 
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FIGURE  I.  Principle  of  AG  focusing.  A focusing-defocusing  or  defociising-focusing  lens  pair  gives 
rise  to  net  focusing.  For  equally  strong  len.ses  the  inter-lcns  distance  should  be  less  than  twice  the 
magnitude  of  the  lens  focal  length.  The  combination  ‘Focusing  (F)-Drift  (0~“/,ero")-Defocusing  (D)- 
Drift  (0~“zero”)‘  is  commonly  referred  to  as  a FODO  cell. 


in  Fig,  2(a),  the  pulse.s  are  applied  alternately  ( fZ  / ||  ) to  the  two  pairs  of  diagonally 
opposite  quadrupole  rods  with  a pha,se  delay  of  T /2,  The  resulting  voltage  difference 
between  adjacent  rods  then  varies  as  shown  in  Fig,  2(c),  The  time-dependent  electric 
potential  near  the  axis  approximates  a perfect  quadrupole 


^x{l)  - <^i{l  - T /2)  - y'^ 

o ..2  ’ 


(1) 


and  gives  rise  to  pseudo-harmonic  radial  confinement.  The  corresponding  radial  oscil- 
lation frequency  (Op  can  be  calculated  from  the  one  period  transfer  matrix  [3,  8],  In 
Eq,  (1),  the  time  elapsed  between  subsequent  focusing  and  dcfocusing  pulses  is  7/2 
and  all  pulses  have  the  same  duration  and  amplitude.  For  an  energetic  ion  beam  in  a 
magnetic  AG  storage  ring,  this  corresponds  to  propagation  through  a regular  so-called 


FIGURE  2.  (a)  The  linear  Paul  trap.  Time-varying  voltages  (t>^(f)  and  ij),(f  -T/2)  are  applied  to  adjacent 
rods  as  shown.  To  obtain  axial  confinement  of  ions,  each  electrode  is  sectioned  into  three  and  an  additional 
dc  voltage  Ude 's  applied  to  the  eight  end-pieces  with  respect  to  the  center-pieces,  (b)  The  pulsed,  time- 
varying  waveform  it>i(r),  (c)  Example  of  the  measured  time-varying  voltage  differences  between  adjacent 
quadrupole  rods  when  the  dual  channel  pulse  amplifier  applies  voltage  pulses  alternately  to  the  two  pairs 
of  diagonally  opposite  quadrupole  rods;  cokf  = 2n  X I MHz  and  pulse  duration  1 00  ns  (i,c,,  T = 0, 1 ), 
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FODO-lattice^.  The  parameters  can,  however,  be  varied  quite  arbitrarily  with  respect  to 
duration  and  delay  and  still  give  rise  to  radial  confinement.  The  pulses  are  produced  by 
a dual-channel  pulse  amplifier  with  inputs  from  two  pulse  generators  triggered  at  the 
desired  repetition  frequency  by  a function  generator.  Each  channel  of  the  pulse  amplifier 
is  rated  to  a pulse  amplitude  of  200  V at  (Orf  = 2jc  x 1 MHz  with  a minimum  pulse 
duration  < 100  ns.  An  axial  confinement  force  is  achieved  by  sectioning  each  of  the 
quadrupole  rods  into  three  as  shown  in  Fig.  2(a)  and  applying  an  additional  dc  voltage 
f/dc  to  the  eight  end-pieces.  This  gives  rise  to  a static,  approximately  harmonic  potential 
*Fj(x,y,z)  = y\ZQ^lJdc[^  — {}?  +y^)/2],  where  T|  is  a geometric  factor.  In  the  current  ex- 
periments, zo  = 5.4  mm  corresponding  to  T|  = 0.248  and  the  harmonic  approximation  is 
valid  up  to  a few  millimeters  from  the  trap  center.  The  single  particle  axial  oscillation 
frequency  is  cOj  = where  a = — 4r|efydc/mzQC0^.  We  note  that  the  axially 

confining  force  arising  from  the  static  potential  T'i  is  accompanied  by  a radial,  defocus- 
ing  force.  Denoting  the  radial  distance  to  the  trap  axis  by  p the  radial  component  of  this 
force  is  = m(co^/2)p. 


Cold  Plasma  in  a Harmonic  Trap 


The  linear  Paul  trap  can  be  described  as  a harmonic  trap  with  an  effective  potential  of 
the  form  (using  cylindrical  coordinates) 


where 


co^  - - (O; 


(2) 


(3) 


Here,  {0<j,  accounts  for  a possible  angular  rotation  frequency  of  the  trapped  plasma  about 
the  z-axis  and  this  defocusing  term  usually  equals  zero  for  a linear  Paul  trap.  We  note 
that  the  meaning  of  P is  nothing  but  the  (resulting)  radial  to  axial  confinement  force 
ratio. 

Neglecting  variations  on  an  inter-particle  scale,  the  density  of  a low  temperature 
plasma  (cold  fluid)  in  a harmonic  trap  is  constant  and  given  by 


no  ■ 


£o 


260^(0^ 


(p+i)= 


2eom , 


U0„ 


■to: 


(4) 


independent  of  (0^.  Furthermore,  the  shape  of  the  cold  plasma  can  be  shown  to  be  a 
spheroid  [9]  and  is  hence  uniquely  described  through  the  length  of  the  envelope  semi- 
axes Z (along  the  trap  axis)  and  R (in  radial  direction).  The  spheroid  aspect  ratio  a = Z/R 
and  the  p-parameter  in  Eq.  (3)  is  related  [9]  as  shown  in  Fig.  3(b). 


^ The  current  setup  is  restricted  to  emulate  the  quadrupole  focusing/defocusing  forces  of  a storage  ring. 
For  more  involved  experiments  on  crystalline  beam  emulation,  one  could  look  into  the  possibilities  of 
implementing  shear  from  bending  by  perturbing  the  crystals  sideways  with  a suitable,  periodic,  transverse 
electric  force. 
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FIGURE  3.  (a)  Example  of  a Coulomb  crj'stal  with  a shape  and  density  described  well  by  cold  fluid 
theory  with  to,),  = 0.  The  ions  are  arranged  in  closed  spheroidal  shells,  (b)  The  relation  between  a and  (3 
as  given  by  cold  fluid  theory  (line),  p for  a string-of-discs  state  crystal  (o)  calculated  from  experimental 
parameters  using  Eq,  (3)  with  0)(,  =0.31  x u>,(. 


The  shape  of  (non-rotating)  Coulomb  crystals  in  a linear  Paul  trap  has  previously  been 
shown  to  described  well  by  cold  fluid  theory  [10].  Figure  3(a)  shows  an  image  of  such  a 
crystal  in  the  pulse-excited  trap. 

Bulk  plasma  o.scillations  in  the  cold  fluid  occur  at  the  plasma  frequency  cOp  = 
\Zc^/io/eo»t  which  also  equals  the  maximum  bulk  mode  frequency  in  the  case  of  a crys- 
tal [1  Ij.  From  Eq.  (4)  we  get  co^  = 2((i-|-  l)o)?  = 2(c0p  — co^).  To  avoid  a 1:2  (parametric 
resonance)  relationship  between  frequency  of  a bulk  mode  and  the  rf  driving  field  we 
have  the  criterion  2(0p  < cOrf.  This  is  equivalent  to  the  tune  condition  for  a crystalline  ion 
beam  as  predicted  by  theory  [12]. 


RESULTS  FROM  EXPERIMENTS 

In  a recent  series  of  experiments,  we  investigated  the  stability  properties  of  laser  cooled 
^‘*Mg"'‘  Coulomb  crystals  when  approaching  the  2(0p  = tOrf  resonance  condition  [13].  In 
this  contribution,  we  shall  only  briefly  summarize  and  discuss  some  of  the  results  from 
this  work.  Notably,  it  was  found  that  trapped  crystals  (initially  in  stable,  closed  shell 
configuration  comparable  to  the  one  .shown  in  Fig.  3)  become  unstable  and  blow  up  to  a 
gaseous  cloud  state  when  the  plasma  frequency  is  increased  to  hit  the  parametric  reso- 
nance. The  plasma  frequency  was  increased  by  increasing  (Op  through  the  rf  waveform 
amplitude  V^f,  while  all  other  parameters  were  kept  constant.  The  observed  stability  limit 
was  in  excellent  agreement  with  the  one  predicted  by  theory  and  seemed  quite  insensi- 
tive to  the  axial  confinement  and  the  number  of  particles  (~50-1000).  Surprisingly,  the 
plasma  was,  however,  observed  to  settle  into  a different  and  hitherto  unobserved  ordered 
structure  when  (Op  was  increa.sed  slightly  above  the  value  of  instability  (see  Fig.  4).  This 
final  state  has  lower  density  and  different  spheroidal  envelope  than  the  initial  closed 
shell  structure  observed  just  before  hitting  the  point  of  instability,  and  appears  as  an  ax- 
ially ordered  string  of  discs.  By  averaging  over  more  frames  we  have,  furthermore,  seen 
indications  of  radial  ordering  within  the  discs  (see  Fig.  5). 

As  mentioned  above,  the  initial  closed  shell  structure  shown  in  Fig.  4 is  described 
well  by  cold  fluid  theory  with  the  plasma  in  a non-rotating  state.  This  cannot  be  the 
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FIGURE  4.  Upper  part:  Frame  sequence  (1-12)  showing  the  closed  shell— *string-of-discs  structural 
transition  for  about  370  ions  when  passing  the  parametric  resonance  2t0p  = cOrf',  the  sequence  was  acquired 
with  5 frames/s  and  the  trap  was  operated  at  Ctttf  = 2ti  x 700  kHz  with  t = 0. 14.  Lower  part;  Close-up 
images  of  the  initial  and  final  state  for  the  structural  transition. 


case  for  the  final  state  as  we  only  increase  slightly  to  induce  the  structural  transition 
so  that  (Op  enters  with  approximately  the  same  amount  in  the  expressions  for  density 
and  shape.  Moreover,  if  one,  anyway,  considers  a slight  increase  in  COp,  this  implies 
an  increase  both  in  density  and  aspect  ratio,  in  contrast  to  the  observations  of  the  final 
state.  However,  the  observed  quantities  might  be  explained  from  the  cold  fluid  equations 
by  assuming  the  plasma  to  be  in  a rotating  state  with  (O^  7^  0 [see,  e.g„  Eq.  (4)].  This 
conjecture  is  supported  by  measurements  performed  on  the  string-of-discs  state  where 
CDj  was  varied  and  the  corresponding  aspect  ratio  was  recorded  [13],  By  assuming  the 
plasma  to  be  rotating  at  the  frequency  0),j,  = 0.31  x (Orf,  we  obtain  the  expected  relation 
between  P and  a as  given  by  cold  fluid  theory  [see  Fig.  3(b)].  Otherwise  (i.e.,  putting 
o),(,  = 0),  P calculated  from  Eq.  (3)  would  be  an  order  of  magnitude  off.  A mechanism 
leading  to  rotation  and  which,  furthermore,  yields  the  assumed  frequency  is  not  clear  at 
present. 


FIGURE  5.  String-of-discs  structure  showing  radial  ordering.  The  image  was  obtained  by  averaging 
over  5 frames  (=1  s). 
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CONCLUSION  AND  DISCUSSION 


In  conclusion,  we  have  investigated  the  stability  of  Coulomb  crystals  in  a pulse-excited 
linear  Paul  trap  having  the  same  transverse  dynamics  as  a storage  ring  with  AG  magnetic 
confinement.  The  observations  agree  well  with  current  theory  on  crystalline  plasmas  and 
beams.  Notably,  it  seems  that  the  criterion  2(Hp  < (Orf  should  be  taken  very  seriously. 
However,  at  the  stability  limit  we  observed  cry.stals  to  survive  by  rearranging  in  axially 
ordered  structures.  The  observation  was  shown  to  be  consistent  with  a string  of  rotating 
discs.  We  recently  learned  about  some  early  work  by  Bliimel  et  al.  [14]  where  similar 
“layered  structures”  were  found  in  connection  with  molecular  dynamics  simulations  of 
up  to  500  laser  cooled  ions  in  a conventional  Paul  trap.  Why  these  simulations  give  ri.se  to 
layered  structures  instead  of  ending  up  in  a closed  shell  structure  as  found  in  equivalent 
work  (see,  e.g.,  [15]  and  references  therein)  is  not  clear  at  present.  The  possible  and 
likely  connection  to  our  work  will  a subject  of  future  studies. 
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